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Analytical Models for 
Battlespace Information Operations 

(BAT-IO) 

Donald P. Gaver 
Patricia A. Jacobs 

Parti 

Some Nonlinear Dynamic Models of Strike 
Response to Region Defense 

EXECUTIVE SUMMARY 

The coordination of information acquisition and interpretation to direct force 
application is increasingly recognized as a crucial military systems design and 
investment issue. This paper illustrates tradeoffs between Blue/own regional 
Attacker sensor and shooter capabilities: it studies a deep strike or SCUD- 
himting scenario in a low-resolution, aggregated manner using an analytical 
state-space approach that recognizes gross aggregated regional Defender (Red), 
and regional Attacker (Blue), system capabilities and limitations. Emphasis is 
accordingly placed on explicitly modeling the availability and utilization of 
information to a striking Attacker, as it becomes available from a realistically 
finite sensor and C2 capability. The (imperfect) information on opposition units, 
the Defenders, that are candidates for prosecution by the Attackers is passed to 
the finite, hence saturable (here missile-firing) Attacking force, the shooters, that 
then responds by prosecuting those units. 



The models specifically recognize that regional Defenders will not be detected 
immediately, nor recognized perfectly, nor are Defender shots (e.g. SCUD 
launcher) fired perfectly, or immediately. Furthermore, attempts to effectively 
target are also realistically modeled as afflicted by imperfect Attacker battle damage 
assessment (BDA), an incapability that, if pronoimced, will non-linearly saturate 
shooters, increase their response times, and hence reduce targeting effectiveness 
and efficient ammunition expenditures. Such models can allow for adaptation by 
both attackers and defenders to recent fortunes: if Defender presence and activity 
is effectively countered by Attackers, then the former may tend to be deterred or 
withdraw; if not, the Defenders are motivated to press their apparent advantage. 
Sharp, threshold-like, responses can follow from the possibly multi-stable 
d)mamics. This behavior will be explored in a second report (Part II). 

The present models are mainly deterministic or pseudo-stochastic in that they 
represent the non-linear effect of stochastic saturation approximately, but 
adequately. However, they can straightforwardly be "made stochastic", 
especially Markovian, and so realized using Monte Carlo simulations. Computer 
programs exist to provide numerical results; some are given. A simple one- 
dimensional stochastic (Markov birth-death) model is given as an appendix to 
Part n. This model can be shown analytically and numerically to exhibit 
"stochastic bi-stability" properties that under certain circumstances (parameter 
combinations) lead to bimodal steady-state distributions. Such a tendency will 
occur also in more detailed,, but less analytically tractable models. 

There are many problera elements that have been initially and purposefully 
ignored. They will be addressed in later work. For instance, the effects of 
different target types, false targets, and decoys must be added (some "decoys" 
are in effect present, in the form of killed Defenders, not so recognized, that are 
mistakenly re-targeted). The effect of different principles for Attacker target 
prioritization under uncertainty, i.e. dynamic scheduling, requires systematic 
attention. In the present models Attackers are invulnerable to attack; this is not 
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always realistic, and can be changed to a duel-like scenario involving 
suppression of enemy (Defender) air defense (SEAD); a paper on this topic is in 
progress. In the current paper Attackers employ generic missiles only, but the 
use of (vulnerable, manned) Attack aircraft can similarly be modeled, as can 
combinations of Attack aircraft. Naval gunfire, and missiles, recognizing the 
coordination difficulties. Emplo)maent of cued reconnaissance aircraft, possibly 
UAVs, can likewise be represented quantitatively as state-space components. In 
addition, refinements that more faithfully represent spatial and perhaps other 
environmental constraints can be incorporated, as can details of conununications 
assets and message-handling protocols in use by both Attackers and Defenders. 

The present papers describe some of the possibilities for insights inherent in 
an enhanced state-space approach. As pointed out, many elaborations are 
possible. The objective is to recognize only that detail in the (preliminary) models 
that is sufficient to hint at payoff from adding suitable assets and strategies at 
appropriate points in the gntire system. Finer detail and resolution is left to 
others to include, and possibly profit by. More elaborate and high-resolution 
models within such tools as NSS (METRON), and JWARS eventually can focus 
with greater intensity on some of the issues raised here. 

In general we belie '/e that this report is in accord with many of the views and 
suggestions of Hachinski (1996), and also of Dockery and Woodcock (1993), and 
others. Those two publications contain many references, some to previous work 
on related topics. 
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Some Nonlinear Dynamic Models of Strike 
Response to Region Defense 
(BAT-IO) 



Donald P. Gaver 
Patricia A. Jacobs 

1. Problem Formulation: Scenarios and Analytical Strategy 

Consider this generic scenario: Defenders (Red units) enter a region, in 
which they assemble, and that they wish to occupy and defend. Their purpose is 
to prepare to oppose the friendly/own (Blue Attacker) assets that, for example, 
may be intending to carry out an amphibious landing. Alternatively, the Red 
Defenders may simply intend to move about H and occasionally shoot harassing 
missiles in the fashion of TEL/SCUD systems at the time of the Gulf War; in this 
case the Blue Attacker facilities may be of lesser and different capability; see the 
OR/MS thesis of Munson (1996). Counter-fire could, in futixre, also at some stage 
come from a Navy force located near, but offshore from . 

The region TV is assumed to be imder Blue surveillance; furthermore, each 
Red shot occasion is an opportunity for detection and response by a Blue 
Attacker. Response to Reds means that the Attacker shoots at the Defender after 
the latter has reached the "head of the line" of a queue or target list of others who 
have come to Attacker attention. Target priorities that order targets in urgency 
classes in the shooting queue are not yet considered here; the effect of such must 
be addressed subsequently. The priorities should be time and experience related, 
and established dynamically to recognize the uncertainty of target list identities. 
Also, all targets in the present model have the same value; this must be subject to 
change. The present model does account explicitly for possible shooter saturation 



or overload; the same is true for the overall Blue sensor force. And the analysis is 
purposefully made time-dependent, although simple steady-state results are of 
some interest to gain insight as to particular combinations of parameters. 

Suppose further that Attacker (Blue) probability of (Red) kill is pK) ideally, if a 
fired-upon Defender is killed, he/she is known to be no longer a threat; while, if 
missed, he/she returns to the pool of Defenders that are candidates for more 
"service", i.e. retargeting. Note that, more realistically, pK might depend on the 
load or backlog of Defender-shooters revealed by Blue's sensors; the sensor, and 
Blue shooters' opportunities are enhanced by the Reds' recent shooting activities. 
The kill probability, px, could plausibly eventually decrease with increases in 
such load because of delays in prosecuting targets: such delays give opportunity 
for the target to move or evade. The initial model treats pK as a constant 
(implicitly depending on only the average range of defender weapon fire that 
reaches the region of attacker location). This feature is a strong candidate for 
refinement in several ways. 

Our models specifically recognize that after a target is fired upon there is a 
battle-damage assessment (BDA) step taken by the Attacker force. Realistically, 
BDA is imperfect. In the present models BDA success is explicitly represented by 
conditional probabilities that reflect the chance of BDA error. One operating 
characteristic of BDA is the conditional probability that the Blue C4ISR system 
concludes that the target is alive, given that the Blue Attacker's response actually 
killed the Red target; denote this by Cda- Also there is a non-zero probability that 
the target is perceived by the Blue BDA system to be dead, given that the attack 
shot actually missed, so the target is actually alive-, denote this error probability 
by Cad- Note that if Cda > 0, as is realistic, an additional — and potentially 
unnecessary — load is gradually imposed on the Blue shooters, and this slows 
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Blue's effective response to Red presence and opposition potential. Of course, 
rapid invasion of the region by Defenders should give rise to many detections of 
Red imits by Blue, which initially produces many opportvmities but, if continued, 
subsequently loads up and may saturate the Blue Attack capability, reducing its 
responsiveness. Consequently it is important to adapt the Attacker force size and 
capabilities to adequately and cost-effectively control a perceived Defense 
intensity in the region. An adaptive control feature that regulates and allocates 
scarce attack assets has been added to our model sequence in Part II. 

A further realistic feature of the present models is the possibility that Red 
units under surveillance may be lost, and hence returned to undetected status. 
Such a loss can be caused by a deliberate Red policy of occasionally leaving 
region IB, or hiding within it temporarily in order to shake off Blue pursuit and 
prosecution. Additionally, losses from surveillance can also occur inadvertently 
because of system inability to maintain effective contact, perhaps as a result of 
terrain properties, but also target action. The present models do not include the 
highly realistic presence of false targets, either of natural origin or deliberately 
introduced by Red as decoys. See Figure 1.1. 

The purpose of this report is to provide first-step tools with which to 
investigate and describe certain tradeoffs and opportunities available in the sort 
of situation described. We begin by analyzing what may be the simplest 
reasonable analytical model(s) that suggest themselves. The models are here first 
exercised in a deterministic, approximately expected-value mode so as to quickly 
and efficiently explore for sensitivities or the "knees in the curves" beloved by 
analysts. In some cases quite simple but suggestive analytical (actually algebraic) 
solutions are possible; more comprehensive time-dependence effects are 
available by package computer programs that solve non-linear differential 
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equations, e.g. MATLAB, or MATHEMATIC A. Stochastic versions of some of the 
models have been formulated, and run in VISUAL BASIC; a simple example is 
available from the authors; it runs under Windows95. Such simulations promise 



Region Defense and Attack 




O Undetected Red Units 

® Detected, Alive (possibly previously targeted) 

13 Detected, Dead (BDA classified as alive; targeted) 

Figure 1.1 



4 



to furnish preliminary rough but economical understanding of the rather 
complex, uncertain, and adaptive phenomena under study. The latter can guide 
further higher-resolution work, if desired. Computer run times are quite short in 
the models proposed, so exploration of alternatives is expedited. 

The modeling strategy espoused is that of first exploring the conflict in an 
aggregated low-resolution broad-brush fashion using simple state-space models. 
This can be done quickly at minimal cost in the scarcest of analyst resources: 
time. Suggestions from such aggregated models can then be investigated at 
higher fidelity and at suitably detailed level using more elaborate models, such 
as perhaps NSS, GCAMS, or the NFS exploratory campaign model JWAEP, and, 
eventually, JWARS. We think of the present low-resolution models as tentative 
and subject to modification (often quickly and easily), both in response to the 
judgments and interests of analysts, and as a consequence of comparison with 
high-fidelity/resolution model runs; ideally, comparisons might even be made 
with exercise or true combat data. We believe that such interplay between model 
levels is a healthy and profitable way to carry out military modeling in support 
of analysis and understanding. Ultimately, models are formulated so as to enrich 
the vision and insights of analysts, at best providing insights at relatively low 
cost in time and money. 
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2. Model I: Surveillance Target Classification (Delayed Saturable BDA) 

Consider this model: Defending Red units enter region at a specified rate. 
They are the subject of a prescribed level of surveillance /reconnaissance, are 
identified (possibly incorrectly) as alive, and are eventually targeted by members 
of a group of Attacking Blue shooters. As a result, some are killed and remain in 
the region imclassified until revisited by the surveillance system; at this point the 
error-prone classification service occurs, and a dead Red may be classified as 
alive and is ultimately (and wastefully) retargeted. Live Red vmits in the region 
are similarly found by the sensor system and targeted, perhaps after several 
previous attempts have been made to shoot them. If a Defender is a TEL, and if it 
fires a (SCUD) missile, it is presumed to be quickly detected (although its track 
may be lost before prosecution). When an unclassified but dead Red unit in the 
region is finally classified as dead that unit is in future considered to be non- 
existent; this may take considerable time, during which it is effectively a decoy. 

Red units may be lost to Blue's surveillance because of a variety of Red 
actions: for example, a TEL may take refuge in a prepared hideout, perhaps 
under a bridge. Such losses to detection and surveillance are assumed in the 
models to occur at a (user-specified) constant rate. Units, once lost, are subject to 
reacquisition at the same rate as newly arrived units. Units in hiding presumably 
cannot fire weapons (laimch SCUDs). 

The setup describes a certain kind of (imperfect) delayed information 
processing concerning the perceived state of both newly-arrived and previously 
engaged live targets, but also the state of targeted Defenders that have been killed 
but not yet classified. 

The above setup represents a system that includes a form of effectively 
delayed BDA; it is one that depends on the general rate of surveillance of the 
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region of interest. A subsequent model addresses the implications of a 
classification system tied more closely to the shooters. 

Parameters of Model I 

A(f) = Defender arrival rate into region . The assumption of a constant 
arrival rate actually permits an explicit steady-state solution to be 
found, for what that is worth. One can perhaps string such together 
to represent various stages of arrival, during which the system 
reaches quasi-stationarity. The model also represents a situation in 
which an initial number of red Defenders is present in R, with no 
reinforcements. 

ji = Service rate of an attacker (1//X = mean time for an Attacker to track, 
shoot, flight time of missile). 

PK = Probability an attacked target is killed. This parameter is presently 
taken to be range- dependent only on average. 

V = Track-loss rate (1/ v = mean of a holding time of Defender in track). 

Once a unit track is lost, it is no longer a viable target for Blue 
service. 

a = Defender (e.g. TEL = SCUD launcher) shoot rate {1/ a = mean time 
between shots by a single Defender). This parameter is irrelevant 
imder certain circumstances. 

7 = Rate Defenders leave region (1/7= mean holding time of Defender 

unit in region), or hide within region. 

^ = Rate undetected /tracked Defenders are acquired by sensors /C2 

system. 

SD = Number of Attacker-shooters; these are viewed as individual servers 
that engage Defender targets for a given time. 

Taa = Probability that a target that is alive is classified as being alive by the 
sensors /C2 system. 1 - = Tad = P(target classified or perceived as 

dead I alive). 

Vda = Probability that a dead target is misclassified as being alive. 

Any or all of the above parameters can be dependent upon time, at the discretion 
of the analyst. We maintain constants for very preliminary steady-state analysis. 
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State Variables for Model I 



The following state variables are needed to describe the present dynamical 
system: ^ 

Ru{t) = Number of undetected live (hence potentially active and threatening) 
Red Defenders present in region at time t 

Rd{t) = Number of detected live Red Defenders present at time t. These are 
on the Blue Attacker-shooter's target list, and will be engaged rmless 
lost by the sensor system (they may go into hiding, or even leave the 
region covered by the surveillance, e.g. JSTARS). 

D(f) = Number of detected and perceived to be alive, hence potentially Blue- 
Attacker-targeted, but actually dead Red Defenders at t. These are 
present because Blue battle damage assessment (BDA) is realistically 
imperfect. 

M{t) = Number of dead Defenders in the region that are not yet classified. 

Classification is done by the Blue surveillance/reconnaissance 
system. 



For the deterministic modeling of saturable service by the Attack shooters 
(presumably missiles in the present context, although Attack aircraft can also be 
represented) we make use of the following approximation, cf. Filipiak (1988), 
Agnew (1976), and also Rider (1967), for the rate of processing by the Blue 
shooters, 

\Rj{t) + D(t)]iiso 
l+[R^(t)+D(()] ■ 
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State Transition Equations for Model I 



dRuit) 

dt 



= ^ + 

Arrival 

rate 



Hc(t)\l-PK] ■ 

Rd{t) + D{t) 

V ' 

Blue service rate of live Defenders 
that results in failure /missed shot 



vRd(0 

Rate of 
track loss 



p'aa^ui^) 

^ V mmm y V 

Defenders leaving region Surveillance /recon. 
and detections caused by system detection and 
Red activity (a) correct classification 

(as alive) 



(2.1a) 



M - 25^) - (2.1b) 

Surv. /Recon. Detection Rate of leaving ' ^ 

detection and by Red region and track Rate of Blue shooting 
classification action loss at live Defenders 

as alive 



Rate of Surv, /Recon. 

detection and 
incorrect classification 
of dead tairgets 



— — — ^s (0 

Rd{t) + D{t) 

Rate of shooter service 
of dead (misclassif. 
as live) Defenders 



(2.1c) 






^dd 



Dead 

classif. 

correctly 

(permanently) 



+ 



g. )M(f) 

Dead 

classif. 

as 

live 






— — — 

Rd(t) + D{t) 

V V ^ 

Service of dead 
but perceived 
as live Defenders 



Ut) 

Rd(f)+D(() 





^s{i)pK 



Service rate of live 
Defenders that results 
in success/kill 



(2.1d) 



Long-Run Solutions 

In order to investigate the possible long-run or steady-state behavior of the 
various state variables, and to compute answers to some interesting operational 
questions, look at solutions to the equations obtained by setting the derivatives 
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equal to zero when X{t) = A, a constant; these are fixed points to which, under 
certain conditions, the system will eventually settle. 



Kd+L> 



0 = + ctR„ - (7 + v)R^ - tt^Hs 

i<d+L> 



0 = 4^d«M- 



D 



Rd +D 



Hs 



Solving for M in (2.2c) and substituting the result into (2.2d) results in 

D = R/f-pK=cRi. 

^dd 

Multiplying (2.2b) by (1 - px) and adding to (2.2a) results in 

0 = A + ( - 7(1 - ))^d - Ru {(XPK + 7 + ^raapK ) • 



Thus, 



•Ku = ^1 + KlR-d 



where 



(XPK+^aaPK+r 

^ ^ V7?K-7(1 -Pk) 

^ CCpK + ^aaPK+r' 

Special Case: 7= 0 (rate at which alive Defenders leave the region is 0) 
If 7= 0, then adding equations (2.2a) - (2.2b) results in 



(2.2a) 

(2.2b) 

(2.2c) 

(2.2d) 

(2.3) 

(2.4) 

(2.5) 

(2.6a) 

(2.6b) 

(2.7) 
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(2.8) 



0 = A- 



1 



1 + c 






1 + (1 + c)k/'= 



where c is defined in (2.3). 



Thus, 



0 = A + A(1 + c)Rd - IJSoRdPK • 

Hence, the condition for the system to be stable is 



(2.9) 



A 



f^DPK 



1 + PK 






\ 



^dd 



< 1 . 



( 2 . 10 ) 



The steady-state number of active Defenders in the region when y = 0 and 
(2.10) holds is 

Ru+Rd=Ki + {K2 + l)Rd 

X 



A ((ajfv)+^ 



^^DPK 



( 2 . 11 ) 



(XPK + ^aaPK + 



aa 



1 -- 



f^DPK 






Discussion 

Examination of the final (right-most) term in (2.11) shows that, if the BDA 
error probability rda is at all substantial, the queue of waiting targets (many 
actually dead) skyrockets non-linearly. This is the result of wasteful retargeting. 
It is therefore a prime technical objective of the Attacker to reduce rda'/ otherwise 
it will be necessary for Attacker/Blue to shoot faster (or with greater effect), 
wastefully using up its ammimition inventory. The same general effect can also 
occur if the Defender region departure or concealment rate, y, is positive but 
relatively small compared to Defender's action rate, a. If this is not so the current 
model would permit Defenders to, somewhat futilely, enter the region and leave 
without shooting. This latter tactic actually might be feasible if only to escape 
possible detection, or, if suspected to be detected, to escape before actual 
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engagement. Or, simply to tie up Attacker forces that might otherwise be used 
elsewhere. 



General Case: y> 0 

If 7> 0, then Rd satisfies the quadratic equation 

0 = (A - >fCi ) + r 4(A - )(! + c) - /isdPk - r(K2 + 1)] - /(Kj + 1)(1 + 

Define 

<j = [(1 - )Ki)(l + c) - USCPK - y{K 2 + l)f + 4(1 - }«i)r(i <2 + 1)(1 + c) 



It can be shown that 



( 2 . 12 ) 



r- 

< (1 + c )- 


+ 7(^2 + 1 ) 


1 4 ( 1 < C 2 + 1)(1 + c ) 


\ / 


L J 


X-jK^ 



(2.13) 



{cc + ^aa)PK + r 



(2.14) 



and that the quadratic equation for Rd, (2.12), will have only one positive root, 
and this provides the long-nm value of Rd- The Ru comes from (2.5), and that of 
D from (2.3). 
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3. Model II: Simplified Surveillance with Immediate BDA 

We introduce next a model that parsimoniously represents salient features of 
the posited long-range Blue response to the Red Defense of ^ assembly 
situation. It is also simple enough to permit explicit algebraic solution for long- 
run or steady-state behavior, if such exists. The latter formulas allow very 
convenient automated exploration by such devices as ANTS (Active Nonlinear 
Tests of Complex Simulation Models), developed by J.H. Miller of Carnegie- 
Mellon University and the Sante Fe Institute. 

In this model there is assumed to be additional ability for the shooter server 
to conduct BDA immediately after firing. 

Additional Parameters of Model II 

Caa = Probability that a Defender that is alive is classified as alive 
immediately after it has been fired upon. I - Caa = Cad = P(Defender 
classified or perceived as dead I alive). 

Cda = Probability that a dead Defender is misclassified as being alive 
immediately after being fired upon. 1 - Cda = Cdd = P(Defender 
classified or perceived as dead I dead). 

The state variables considered are these; 

State Variables for Model II 

We need the following variables to describe system evolution; 

Ru{t) = Number of undetected live (hence potentially active and threatening) 
Red Defenders present in region at time t. 

Rd(t) = Number of detected live Red Defenders present at time t. These are 
on the Blue Attacker-shooter's target Ust, and wiU be engaged unless 
lost by the sensor system (some may go into hiding, or even leave 
the region covered by Blue surveillance, e.g. JSTARS). 

D(f) = Number of detected and perceived to be alive, hence potentially Blue- 
targeted, but actually dead Red Defenders at f. These are present 
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because Blue battle damage assessment (BDA) is realistically 
imperfect. 

In this model Blue Attackers ignore the Red Defenders that are dead and 
once-perceived to be dead; if a dead target is ever so classified, it is thereafter 
omitted from consideration. In this model, classification (BDA) is modeled as 
carried out soon (immediately) after a Blue Attacker engages /shoots at a Red 
Defender. 



State Transition Equations for Model II 

The following are transition equations for evolution of the system state. 



dRjM ^ ; l (0 

dt ‘-V— 

Arrivals 
into region 

— — . : ' Miss Mis-class 

Blue service / shoot rate (at live Red target) Probability ("dead" back 

to undetected) 



+ Rd{^)v Ru {f)fyaa R^ui^)^ ^uij)y 

Alive lost by New New detections Red Defenders 
surveillance detections caused by leave region 
system (sensors) Red activity 

(shooting) 



New sensor Activity Leave Alive lost by 

detections detections region surveiUance 

system 



MU — H(RM) + D{tyM,Sn)- 

Rrf(f) + D(f) 

y / 

Serv./ shoot rate 



[(l-PicKd + ^l 

Miss & 
misclass 



(3.2a) 



(3.2b) 
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N ^ ^ 

Blue kills an cilive Red, 
misclass. as alive 

' V ' 

"Re-kill" dead classed alive; 
class as dead 



where 



H(Rj(()+D(();m,sd) 



[Rrf(i) + D(t)]/lsp 
l + [Rj(i) + D(()] ■ 



(3.2c) 



Solutions 

In order to investigate the long-run behavior of the various state variables, 
and to compute answers to some interesting operational questions, look at 
solutions to the equations obtained by setting the derivatives equal to zero when 
X{t) = A, a constant; these are the fixed points to which, under certain conditions, 
the system will eventually settle. 

° ^ ^ +D}^i,sj:))+vRd-{a + r)i^ - ^aa^u (3-3a) 

Rd + D 



0 = ocR^ + -{y + v)Rd - ■ + PK^Rd + D;/x,sd) 

Rd + D 

°^J7^[PK‘daMRd +D:fl,SD)-J^CMH(Rj*D:ll,SD) 



where 



<7JC = 1 - PK- 

If the third equation (3.3c) is solved first we find 

^_ PKi'^-Cdd)Rd 

<^dd 



(3.3b) 



(3.3c) 



(3.4) 
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so that 



R. 



■d _ 



+ D 2 ^ PK{'^-<^dd) ' 

^dd 



(3.5) 



To eliminate the nonlinear term from the first two equations (3.3a) - (3.3b) 
multiply the first by {(^K^ad + ptd arid the second by qK^ad and add to obtain 



0 = [X + vRd- {^r^a + a + y)Ru l^K^ad + PK ) + [i^raa + Cc)Ru -{y + v)R^ ]qKCad ■ 
Simplify: 

0 = ^^K^ad + Pk) + i^K - m(^ad)Rd ~ [m^^ad + {^aa + « + 7]Pk]Ru ■ 

Thus, 

^ ^‘iKCad+PK) , in’K-m^^ad) j, 

“ mCad+{^raa+a + r)pK mCad+i^aa + CC + r)pK 

= L + MRa. 

To find Rd, insert (3.9) and (3.4) into (3.3b) and solve for Rd- 

Special Case: 7= 0 (the rate at which alive Defenders leave the region = 0). 
Let 

Cdd 

Assume 7=0. 

Adding equations (3.3a) - (3.3b) results in 

= A[l + (1 + c)Rd] - {JS^ pj^R^ 

= X + Rd{^l + c) - psj^ pi^}. 



(3.6) 



(3.7) 



(3.8) 

(3.9) 



(3.10) 



(3.11) 
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Hence, 



0 = A H 

1 + c 



<lKCad^^D 



{l + c)Rd 
l + {l + c)Rd 



1 

1 + C 



{(jKCad+PK)MSD 



(1 + c)Rd 
l + {l + c)Rd 



= A - /zsd pj^ 



Rd 

l + {l + c)Rd 



= A[l + (1 + c)R^ ] - pi^R^ 



= A + Rrf{A(l + c)-/iSD/?K}- 



(3.12) 



Thus, for the system to be stable 



A 

U^D PK 



1 + PK 



0--Cdd) 
^dd . 



< 1 . 



(3.13) 



In the unrealistic but instructive case in which Defenders never leave the 
region we can thus see directly the effects of BDA misclassification. The steady- 
state number of active Defenders in the region when 7 = 0 is 



Ru+Rd — ^ + (-M - 1 - 1)R^ 



_ ^(^KCgd + PK) . {y + ^aa+Cc) 



+ ■ 



P^DPK 



(3.14) 






1 -- 



P^dPk . 






^da 

^~^da 



We see from the last term in (3.14) that increases in Cda can result in more sizable, 
rapid, and non-linear increases in the total active Defenders, R« + Rd, than do 
increases in Cad) the latter is also influential but only linearly. This behavior 
occurs because Cda > 0 results in unnecessary work by the shooting server and 
can result in system overload. The C 2 server is recognized to be equivalent to a 
deterministic infinite server, which cannot be saturated — a feature that is a 
candidate for change. This model differs only slightly from Model I but it is 
simpler. 
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Comparison of (3.14) with (2.11) in the case in which r^a = Cda shows that the 
alive Red population, Ru + Rd, is smaller for the present immediate BDA model. 
Model n. In this case Ru is smaller for Model II than Model I, while Rd is the same 
for both models. 



General Case: y> 0 

If 7> 0, then Rd satisfies the quadratic equation 

0 = (A - ji) + [(A - )L)(1 + c)-iisd Vk- - r(M + 1)(1 + c)rI (3.15) 

Define 



<7 = [(A - ■)£)(1 + c) - /jsq - y{M + 1)]^ + 4(A - yL)y{M + 1)(1 + c) 



= {X-yLY 



It can be shown that 



(1 + c)- 



A-tL = 



MSdPK+7{^ + ^) 
(A-)L) 






n2 



+ 4 



(M + l)(l + c) 



(3.16) 



A-)L 






> 0 , 



(3.17) 



mCad+{^raa+a + r)pK 
and that the quadratic equation for Rd (3.15) will have one positive root. That 
positive root is given by 



R^ = 



_(A-)L)(l + c)-/iSppK-y(-M + l) + V? 



(3.19) 



2r(M + l)(l + c) 

Note that c is given by (3.10), M by (3.8), q by (3.16). From this formula we can get 
an explicit algebraic expression for Ru from (3.8), and for D from (3.4). 

Numerical Examples 

Two measures of evaluation (MOE) are the rate of Defender attrition 

Rd Rd "I" D Rd 

Pa = „ ^ ■ n , t^ p^dPk = mr—^P^DPK 



Rd + D l+Rd+D 



l + Rd+D 
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and the rate of Defender firing (shooting) 



(p={Ru+Rd) (X. 

Figures 3.1 - 3.11 present results of the two models. Figures 3.1 - 3.10 plot rate 
of Defender firing, (p, versus rate of Defender attrition, pR, for various values of 
a, the firing rate per Defender. The value of each rate is plotted with a symbol for 
each value of a. The model parameters for Figures 3.1 - 3.6 are X = 10, pK = 0.7, 
7= 0.5, v= 20, rda = rad = Cda = ^ad = 0.3. In each figure the rates of firing per 
Defender, a, are 0.1, 0.5, 1.0, 1.5, ..., 4.5, 5. In Figures 3.1 - 3.3, the rate of 
Attacker-shooter service, jiSD = 25. In Figures 3.4 - 3.6, the rate of Attacker- 
shooter service, jusp = 50. In Figures 3.1 and 3.4, the Attacker-sensor acquisition 
rate <^ = 1. In Figures 3.2 and 3.5, the Attacker-sensor acquisition rate (^ = 5. In 
Figures 3.3 and 3.6, the sensor acquisition rate ^ = 10. 

The immediate BDA model always results in a higher rate of Defender 
attrition for the same rate of Defender firing. However, the difference becomes 
negligible for large the rate at which sensors acquire the Defenders. The 
difference is larger when the rate of shooting service, psu = 50, (in Figures 3.4 - 
3.6) than when psD = 25, (in Figures 3.1 - 3.3) for the same values of This 
behavior results from the increased number of xmdetected Defenders, Ru, in the 
delayed BDA model. Note that for small the increase in the rate of Defender 
attrition is more responsive to change in the firing rate per Defender, a; 
Defenders are then more often discovered when they reveal themselves. For 
large surveillance rate, the rate of Defender attrition is almost constant as a 
function of a. Not surprisingly, for small a Defender is more likely to be 
detected just after it fires than before. 

In Figure 3.7, v = 0, | = 1, psD = 10 and ps-Q = 25 with the other parameters the 
same. Note that when psD = 10, the difference between immediate BDA and 
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delayed BDA is negligible. This is because the Attacking-shooting server is 
saturated although the acquisition rate ^ is relatively small. The effect of 
increasing the missile firing rate to piso = 25 is to decrease the rate of Defender 
shooting by decreasing Rd) further, since the shooting server is not now saturated 
there is a larger difference in the rate of Defender attrition between immediate 
and delayed BDA. 

In Figure 3.8, ^ = 1, }^d = 25 and v, the rate at which detected Defenders are 
lost from track, is set equal first to 0 and then to 10. Note the anticipated higher 
rate of Defender attrition for v = 0, and also the accompanying decrease in rate of 
Defender shooting attributable to the decrease in Ru + Rd> the Defender 
population available to fire missiles. 

In Figure 3.9, v = 0, /xsd = 10 and <^, the rate of surveillance /recormaissance, is 
first set equal to 1 , and then raised to 20. There is no apparent difference between 
immediate and delayed BDA for ^ = 20; live targets are reacquired relatively 
quickly. Further, as a increases, the saturation of the shooting server makes the 
rate of Defender attrition equal for ^ = 1 and ^ = 20. The present Blue Attacking 
shooter force can not profit by the increased acquisition capability. 

In Figure 3.10, <^ = 10, v = 20, A = 10, 7 = 0.5, pK = 0.7 and jdSD = 25. One curve, 
with larger rate of Defender attrition, corresponds to rad = ^da = Cad = Cda = 0.3. 
The lower curve corresponds to higher error probabilities, arbitrarily rad = Ua = 
Cad = Cda - 0.5. The less-effective BDA results in a decrease of about 5 in the rate of 
Defender attrition and an increase of about 10 in the rate of Defender shooting 
rate for a = 4.0, 4.5 and 5. The difference between immediate and delayed BDA is 
small because of the comparatively high rate of Blue surveillance/ 
recormaissance, (^ = 10 . 
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In Figure 3.11, A = 10, <§ = 1, = 0-7/ v = 20 and usd = 25. The rate at which 

Defenders leave the region is 7= 0a where a = 1 and 0= 0.1, 0.2, 0.5, 1, 2, ..., 10. If 
0 < 1 then delayed BDA results in a larger Defender shooting rate than 
immediate BDA. As the rate at which Defenders leave the region becomes larger 
than the firing rate per Blue Attacker, the rates of Defender shooting for 
immediate BDA and delayed BDA become essentially the same. The rate of 
Defender attrition is always below the rate of Defender shooting but the rates 
become comparable as the rate of Defender departure from the region becomes 
greater. 



Numerical Example: Nonstationary Results 

In this example there is a maximum number of Red Defenders, M, which 
enter the region at a linearly increasing rate, A(t) = At, for t < Tm where Tm = 



2M 



; note that the total number of arrivals is Asds = M. 

JO 



In Figures 3.12 - 3.15, the total number of Defenders M = 50 and A = 5; thus, 
Tm = is the time at which all Defenders have entered the region. Other 
parameters are as follows: the track loss rate v = 0.5; the detection rate = 0.1; the 
Red activity rate a = 0.5; the Blue service rate jjsd = 12; the probability of kill 
PK = 0.7; the probability a live Red target detected by a Blue sensor is classified as 
live, Taa = 0.8; and the probability Blue classifies an alive Red as dead, Cad = 0.7. 
The Defenders do not leave the area; 7= 0. 

In Figures 3.12 - 3.13, the probability with which Blue correctly classifies a 
dead Red as dead is c^d = 0.2. This rather low figure is reflected in a dramatic 
growth in the backlog of already-dead targets and in a correspondingly drawn 
out campaign. In Figures 3.14 - 3.15, Cdd = 0.8. Figures 3.12 and 3.14 present the 
number of alive Red and dead Red targets (but not yet classified as dead) waiting 
for or being served by the Blue shooters for Model n. When Cdd = 0.2, the number 
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of dead Red targets that are misclassified as alive and are awaiting retargeting is 
much larger than the number of live targets waiting to be served by the Blue 
shooters. The Blue shooters are saturated by the misclassified dead Red targets. 
Figures 3.13 and 3.15 present the accumulated number of Blue shots. Red shots, 
and Reds that are killed for Model II; also presented is the number of Red 
Defenders alive at time t. Comparison of Figures 3.13 and 3.15 indicates that the 
saturation of the Blue shooters by dead Red targets when Cdd = 0.2 decreases 
Blue's ability to prosecute alive Red targets. This impairment results in Red 
Defenders being alive for a far longer period of time and being able to shoot 
many more times. Further, Blue wastes much ammunition retargeting dead Red 
targets. 



22 



4. Model III: Command and Control Delays Explicitly Represented 

The previous models are next extended to reflect realistic delays in 
classification and communication of defender detection to shooters. Although 
delays were modeled before, the very realistic effects of non-linear congestion 
and queuing was not explicitly represented. It may be seen that additional state 
variables are now required to minimally specify the dynamical state of the 
system. As before, 

Ru(t) = Number of undetected live (potentially active) Defenders present in 
the region at time t, 

Ra{t) = Number of newly-acquired live Defenders present at time f, 

Rd{t) = Number of detected live/ functional Defenders present at time t, 

D{t) = Number of detected and misclassified as live, but actually dead. 
Defenders present at time t. 

The Defenders enumerated by Rd{t) and D(f) are viewed as targetable, i.e. 
eligible for Attacker engagement; they are effectively "queued up" for shooter 
service. Additionally, we wish to define and enumerate those Defenders that 
have been engaged and await classification as to damage status (BDA) and 
possible attacker response or service: 

Sfl(f) = Number of live Defenders that have been engaged /shot 

at by Attackers, at t, 

Sd{t) = Number of dead unclassified Defenders present at time t. 

Note that all such variables might be treated as state variables of a multidimen- 
sional birth-death Markov stochastic process. A mathematically explicit 
treatment of such a setup appears cumbersome, but simulation models have been 
written to allow the effects of randomness to be investigated. 
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Parameters of the Model 



T = Maximum number of Defenders allowed in the region 
simultaneously by the Defender decision makers. This parameter is, 
in effect, a control variable. 

X = Arrival rate of Defender (targets) to area. Note that the overall 

arrival rate is here allowed to depend on a hypothetical goal for the 
Defenders, T. Additions are made on the basis of current Defender 
count deficiency from the goal level. 

^ = Rate at which Defenders are "served" by an Attacker-shooter. 

SD = Number of shooters; these are viewed as individual servers that 

engage Defender targets for a given time. 

T] = Rate at which Defenders are served by the Attacker's C2 system, 

viewed as a (saturable) service subsystem. 

su = Number of C2 servers possessed by the Attacker force. 

V = Individual rate at which Defenders (Red targets) are lost from track 

but are still in the region. 

a = Individual rate at which Defenders are active (e.g. shooting), hence 

causing potential damage but also revealing their presence; 
equivalently, 1/ a is the expected time between shots by a Defender. 

^ = Rate of acquisition of Defenders by Attacker surveillance (a sweep- 

width concept). 

y = Rate at which Defenders leave the region; equivalently 1/y is the 

expected time that a Red Defender spends in the region (if it is not 
killed). 

PK = Single-shot probability with which a Blue Attacker kills a Red 

Defender. 

Caa = Probability that an engaged Defender that is still alive is correctly 

classified as being alive; Cad = 1 - Caa is the probability of misclassi- 
fication of a live Defending target as dead. 

Cda = Probability that an engaged Defender that has been killed (is dead) 

is misclassified as being alive; Cdd= 1 - Cda is the probability of 
correctly classifying a dead Defending target as dead. 
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Saturable C2 Service Submodels 

For the deterministic modeling of saturable service by the C2 facility and by 
the Attacker-shooters (presumably missile launchers in the present context) we 
make use of the following approximations, (again Fdipiak (1988)): for the C2 rate 
of service, i.e. of refining and processing information concerning those Defenders 
detected, we use 

Hc(f)HH(K„(0 + S,(f) + Sd(f);T7,sc) 

and for the rate of processing by shooters, 

Hs(f)HH(Rrf(0 + D(f);/t,SD) 

\Rd{t) + D{t)]usD ^ (4.2) 

l + (0 + ^(0] 

Both of these expressions reflect the saturability of the respective service systems. 
It is argued, e.g. for the C2 service system, that the load or backlog on the system 
is the sum of (a) Ra{t), the number of acquired Defenders never before engaged; 
(b) Sa(t), these previously-engaged and alive but unclassified at t; and (c) Sd(t), 
these Defenders previously engaged that are dead but unclassified at f. These 
constitute the entire load seen by the C2 system; that load is processed at 
individual rate 77, so for small load the gross processing rate is proportional to that 
load, while if that load grows the gross service rate saturates at the overall gross 
rate rjsc, as it should. An identical argument holds for the shooter service facility 
behavior. 

We also incorporate a static (time and event independent) weighting/ control 
scheme that prioritizes the server attention to the different classes of targets. 
Again for the C2 service, we allow different C2 service emphasis on, for instance. 
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the newly-detected Defender targets that have not yet been engaged than on the 
previously engaged, but not-yet-damage-classified. The weights can be made to 
adapt to changing backlogs and opportunities, i.e. to represent "emerging 
behavior", but no details are given at present. 



State Transition Equations 

Here are the state transition equations proposed to represent the enhanced 
model. 



dRujt) 

dt 






1- 



Ru{^) + Ra{^) + Rd{^) + ^g{^) 



Defender entry rate to region 
(or emergence from hiding). 
Controlled by Red 



' V ^ 

Rate that unclassified live targets 
are designated by C2 as dead 
(back to undiscovered status) 



+ v[R,{t) + RM + SM) - (y + o + |)R„(f) 

V— — V ' V ' 

Rate of loss of track of Rate of loss: leave region (^); 

Defenders under track at t detection as a consequence of 

Defender action (a); detection of 
Defenders from surveillance action (<^) 



( 4 . 3 ) 



Detection of Red Rate of leaving 

Defenders by their action (a) region (y) or track loss 
and by Blue surveillance ((^) or temporary hiding (v) 

^ARqjt) „ i,s 

u.„R„(t) + »s(S,(() + Sj(t)) ’ 

' V ' 

Rate of C2 system classification 



( 4 . 4 ) 



dRdi^) _ ^ARg{^) ij / . 

' v — ^ 

Rate of C2 system classification of live targets as alive 



U‘) 

R^(0+D(f) 



Hs(t) 



Rate of engagement of 
live, classified Defenders 



{v + y)Rd{t) 

' ^ ^ 

Rate at which track 

lost (v) or region left (y) 
by live classified Defenders 



( 4 . 5 ) 
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‘iSait) 

dt 



Rate at which eligible live Defender targets 
are engaged and missed, becoming 
candidates for classification 



jj- /A 

V ^ . 

Rate of correct classification of 
live unclassified Defenders 



- (v+r)Sfl(0 

' V ' 

Rate of track loss (vl 
or leavmg region (r) 



(4.6) 



rfSd(l) 

dt 



V ^ , 

Rate at which eligible live 
Defenders are engaged and killed 

^ 

(0 (Sfl (0 + (0) 

Rate at which unclassified dead Defenders 
are (correctly) classified as dead 



D(f) 






Hs(0 (4.7) 



Rate at which dead, but 
misclassified as alive. Defenders 
are engaged, rejoining the C2 
classification queue 



dD{t) 

dt 






— ' — ■ V 

Rate at which Defenders that are dead but 
misclassifed as cdive are serviced by C2 and, 
again, misclassified as alive 



D(t) 



U‘)+m) 



Hs(f) 



Rate of engagement of 
Defenders that cire dead, 
but misclassified as alive 



(4.8) 



Numerical Examples 

Table 1 presents steady-state results of Model III with parameters T = 100, 
A = 10, pK = 0.7, V = 0.5, Y = 0.5, a = 1, Cad = Cda = 0-3/ tfs = = 1- Displayed are 

rates of Defender shooting, p$, and rates of Defender attrition pA where 

Ps = a{Ru + Ra+Rd + Sa) 

Note that increasing the acquisition rate ^ from 1 to 10 when rjsu = 10 and 
ps£) = 10 or 20 does not change the Defender shooting and attrition rates by 



27 



much; this behavior is due to saturation of the C2 server. Increasing that capacity, 
as measured by to 30 allows the Defender shooting and attrition rates to 
become more sensitive to changes in ^ and fJSD. 

TABLE 1 

Defender Shooting Rates and Defender Attrition Rates 
for Model with C2 Server 



T = 100, X = 10, PK = 0.7, V = 0.5, 7= 0.5, a = 1, Cad = Cda = 0«3/ '^S = ‘^A = 1 



1 


1 


10 


IJSD 


10 

Defender Rate of 


20 

Defender Rate of 


10 

Defender Rate of 


20 

Defender Rate of 




Shooting 


i Attrition 


Shooting 


Attrition 


Shooting 


Attrition 


Shooting 


Attrition 


psu 

10 


12.7 


2.4 


12.4 


2.6 


12.6 


2.5 


12.3 


2.6 


20 


10.3 


3.8 


9.1 


4.6 


9.8 


4.1 


8.4 


5.0 


30 


9.7 


4.2 


7.6 


5.5 


8.9 


4.7 


5.6 


6.6 



5. Concluding Discussion 

The equations describing the models in this report have been solved 
numerically, beginning with various hypothetical initial conditions, and the 
results are intuitively agreeable. Various functionals of these solutions that are 
operationally meaningful and informative can also be numerically evaluated; e.g. 
the (expected) number of weapons fired by each side up to tie f; this is the 
important MOE inventory expenditure. Cumulative attrition can also be computed, 
as well as the effect of attacks SCUD firings. It is the dependence of such MOEs 
upon system properties, such as surveillance capability, including the adequacy 
of BDA assessment, and attrition capability by defense shooters that the models 
will help to clarify. Tradeoffs can be vmcovered, and the operational value of 
certain proposed system enhancements revealed. We plan to cover this 
exploration phase more extensively in later work. 
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In Part II of this report series we study the effect of adaptive policies that 
regulate both attacker and defender changes in region occupancy as a 
consequence of experience. Some such control phenomenon is dependent on 
perception, which may well be faulty. In subsequent work still more realistic 
issues are included. 

Finally, it is emphasized that this modeling effort is purposefully aggregated 
and broad-brush so that quick turnaround results are possible. The deterministic 
models presented here can be very quickly exercised on modern PCs. The 
simulation models that follow these are similarly quickly run. Such a feature 
means that considerable exploration is possible before the high-resolution 
digging begins. We are convinced that the preliminary spadework described can 
efficiently locate paydirt for later exploitation. 
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RATE OF DEFENDER ATTRITION VS DEFENDER FIRING 
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RATE OF DEFENDER ATTRITION VS DEFENDER FIRING 
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RATE OF DEFENDER ATTRITION VS DEFENDER FIRING 
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RATE OF DEFENDER ATTRITION VS DEFENDER FIRING 
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RATE OF DEFENDER ATTRITION VS DEFENDER FIRING 
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RATE OF DEFENDER ATTRITION VS DEFENDER FIRING 
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RATE OF DEFENDER ATTRITION VS DEFENDER FIRING 
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RATE OF DEFENDER ATTRITION VS DEFENDER FIRING 
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